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ABSTRACT 

Context. To understand the mechanisms that trigger solar flares, we require models describing and quantifying observable responses 
to the original energy release process, since the coronal energy release site itself cannot be resolved with current technical equipment. 
Testing the usefulness of a particular model requires the comparison of its predictions with flare observations. 

Aims. To test the standard flare model (CSHKP-model), we measured the magnetic-flux change rate in five flare events of different 
GOES classes using chromospheric/photospheric observations and compared its progression with observed nonthermal flare emis- 
sion. We calculated the cumulated positive and negative magnetic flux participating in the reconnection process, as well as the total 
reconnection flux. Finally, we investigated the relations between the total reconnection flux, the GOES class of the events, and the 
linear velocity of the flare-associated CMEs. 

Methods. Using high-cadence Hor and TRACE 1600 A image time-series data and MDI/SOHO magnetograms, we measured the 
required observables (newly brightened flare area and magnetic-field strength inside this area). RHESSI and INTEGRAL hard X-ray 
time profiles in nonthermal energy bands were used as observable proxies for the flare-energy release rate. 

Results. We detected strong temporal correlations between the derived magnetic-tlux change rate and the observed nonthermal emis- 
sion of all events. The cumulated positive and negative fluxes, with flux ratios of between 0.64 and 1.35, were almost equivalent to 
each other. Total reconnection fluxes ranged between 1.8 x 10"' Mx for the weakest event (GOES class B9.5) and 15.5 x 10^' Mx 
for the most energetic one (GOES class X17.2). The amount of magnetic flux participating in the reconnection process was higher in 
more energetic events than in weaker ones. Flares with more reconnection flux were associated with faster CMEs. 

Key words. Sun: flares - Sun: corona - Sun: chromosphere - Sun: magnetic fields ~ Sun: activity 



O 1. Introduction the CSHKP model ('Carmichaefl 964": ISturrocklll 9661: iHiravamd 
T-H 11974 : Kopp & Pneuman 1976). According to this, the separat- 
Solar flares are intriguing, intensely-studied phenomena. Much j„g fl^^g jj^bons are flie chromospheric signatures of the energy 
O progress has been made m understandmg the processes that oc- release in the corona. The energy, which has been released at 
L: cur on the Sun when a flare occurs. However, these events re- t^e primary energy-release site, is transported downward into 
^ mam a pardy-unsolved mystery. It is generaUy accepted that the j^e chromosphere along the newly reconnected field lines, by 
[X^ energy fuehng a solar flare is magnetic energy, which is accu- f^^t particles and thermal conduction. When it is deposited at 
^ mulated slowly and stored inside the corona. After its release footpoints of the field lines, which are rooted inside op- 
, by magnetic field reconnection, it is converted into the kinetic ^^^^^^ magnetic polarities, the chromosphere brightens, i.e., the 
energy of fast particles, plasma flows, heat, and MHD shock so-called flare kernels are created, and many of these adjacent 
waves. Unfortunately, our observations remain insensitive to the j^gj-ngj^ f^^^ elongated flare ribbons. The reconnection site is 
primary energy-release process, since the energy release site in ^^^^^^^ ^elow the flare-associated erupting filament/CME, and, 
the reconnecting current sheet cannot be resolved with our cur- gin^e the erupting structure rises continuously to higher coro- 
rent technical equipment. To learn about the properties of the ^^j heights, the reconnection also occurs at successively higher 
energy release site and the processes occurring there, we there- altitudes, as the flare progresses. As a consequence, field lines 
fore depend on exploring the plasma responses to the original ^j^h footpoints rooted further and further away from each other 
energy-release process. are involved in the reconnection process. Thus, the newly bright- 
One of the observable consequences of the energy release ened areas of both ribbons are located further and further apart 
in the corona are the bright, sepai-ating flare ribbons, which and we observe separating flare ribbons, 
can be seen in UV bands and chromospheric spectral lines, 

and are prominent in Ha. They are located on each side of the Since the coronal reconnection site is coupled with the chro- 
magnetic polarity inversion Une and move away from it dur- mosphere by the reconnected field lines, it is obvious that the 
ing the course of the flare. This apparent motion can be ex- coronal magnetic reconnection rate, i.e., the rate at which mag- 
plained by the most widely accepted model of erupting flares, netic field lines are reconnected, is associated with the separat- 

i ng fla re rib bons. Based on the C SHKP model, Forbes & Pries^ 
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tion, they derived a quantitative estimate for the coronal recon- 
nection rate or magnetic flux change rate ip, respectively, that can 
be derived from observable quantities 



da. 



(1) 



where B„ is the photospheric magnetic-field strength component 
perpendicular to the solar surface inside the newly brightened 
area da, which is swept up by the flare ribbons. We denote by 
ip the rate at which magnetic flux from formerly separated mag- 
netic domains is unified. The original CSHKP model is a 2D 
model, since it describes the evolution of a flare in a vertical 
plane. In the third dimension, namely, the direction of the po- 
larity inversion line, translational symmetry is assumed. This as- 
sumption converts the model to a 2|D model. It is likely, how- 
ever, that the extension in the third dimension is not continuous, 
but rather highly fragmented into temporary magnetic islands. 
However, we note that calculating the magnetic flux change rate 
ip does not require the assumption of translational symmetry. 
Therefore, (p is valid in three dimensions. 

To test the CSHKP model, an observable quantity is re- 
quired that can be regarded as a proxy for the energy release 
rate or reconnection rate, respectively, which we are incapable 
of measuring directly. The fast electrons that spiral downward 
along the newly reconnected field lines, generate microwave 
gyrosynchrotron emission, and when they deposit their energy 
at the footpoints of the reconnected field lines, hard X-ray 
(HXR) emission is produced by the nonthermal bremsstrahlung 
of electrons being scattered off ions. Therefore, the observed mi- 
crowave and HXR fluxes act as indicators of the rate of acceler- 
ated electrons, which carry a large fraction of the total energ y 
released during a flare (e.g.. lHudsonlll99itlDennis et alj|2003h . 
Thus, microwave and HXR emission can be used as proxies for 
the flare energy-release rate. If the CSHKP model is applicable, 
the derived magnetic-flux change rate (f should exhibit a similar 
temporal evolution as the observed HXR and microwave flux, 
i.e., the maxima in the (^-profile should coincide approximately 
co-temporally with peaks in the observed nonthermal flare- 
emission. Strong temporal correlations of observed nonthermal 
emission with derived quantities characterizing the magnetic re- 
connection process, such as the magnetic flux change rate, the 
electric field at the reconnecting X-point, or the Poynting flux, 
which is transporte d into the reconnection region, were f ound in 
sever a l cases (e.g . . lAsai et al. 2004: Qiu et al. 2004; Jing^t alj 



20051; llsobe etall 120051: iLee et all 120061: ll^meret alj l2007t 



Miklenic et alj|2007 l)" 



It is known from observations that the energy release in a dy- 
namical fla re is also related to the kinematics of the associated 
CME (e.g.. iMaricic et al.ll2007l : iTemmer et ani2008) . Statistical 
studies indicate that CME parameters, such as the velocity and 
kinetic energy, are correlated with the soft X-ray (SXR) peak 
flux and the time-integrated SXR flux, i.e., the flare fluence 



(e.g..lMoon et alJl2002rrVrsnak et alj|2005h . lOiu & Yurchvshvnl 
reported that a greater amount of total reconnection flux 
is related to higher CME velocities. 

In this paper, we present a detailed analysis of five well- 
observed erupting flare events of different GOES classes. We 
calculated ip to test the model prediction, namely, the co- 
temporal evolution of the derived magnetic-flux change rate and 
observed nonthermal flare emission. In addition, we determined 
the cumulated reconnection-flux profiles for the positive and 
negative magnetic polarity domain as a function of time, the ra- 
tio of these fluxes, as well as the total flux that had to have been 
reconnected until the end of the flare. Finally, we investigated the 



relations between the total reconnection flux, the GOES class of 
the events, and the linear velocity of the flare-associated CMEs. 

The paper is structured as follows. Section |2] contains de- 
scriptions of the used data sets. Sect. [3] describes the methods 
that we have applied. Sect. |4] lists the results for each flare, and 
in Sect.|5]the results are summarized and discussed. 



2. Data and observations 

The five flare events presented in this paper were eruptive (two- 
ribbon) flares, associated with a CME. Table [T] overviews the 
flare/CME events. In each flare, two ribbons are clearly shaped, 
and in some cases also a third ribbon or distant brightenings can 
be seen. The position of a flare on the solar disk served as a selec- 
tion criterion, since the magnetic field component normal to the 
solar surface is required for the analysis (see Sect. [3]). Therefore, 
we chose only events that were located not more than ±40° away 
from the disk center A second selection criterion was the avail- 
ability of high-cadence image time-series data in Ha or UV and 
complete coverage of the impulsive phase in HXRs. 

To analyze the flare events, we used the following data sets 
(see also summary in Table|2]l: 

1. Full-disk line-of-sight mag netograms, pro vided by the 
SOI/MDI instrument (Sche rrer etal] Il995h onboard the 
Solar and Heliospheric Observatory (SOHO). In all events 
except for that on 2003 November 18, five successive mag- 
netograms one-minute apart were averaged together to form 
a single low-noise magnetogram, which was used in further 
analysis. 

2. High-cadence (4 - 60 s) image time-series data in Ha or UV. 
The Ha images were provided by various ground-based ob- 
servatories: Kanzelh ohe Solar Observatory (KSO, Austria; 
'Otruba & Potzi"2003), Big Bear Solar Observatory (BBSO, 
USA; Denker et al. 1999), Hvar Observatory (Croatia; 
lOtrubal l2005h . and Meudon Observatory (France). UV 
images in the 1600 A passband were obtained from 
the Transition R egion and Coronal Explorer (TRACE; 
lHandveta ri999'). 

3. Full-disk, nonthermal HXR-intensity time-profiles from the 
Reuven Ra maty Hig h Ener gy Solar Spectroscopic Imager 
(RHESSI; iLin et alj llOOl) . For one event (2003 Oct 
28), HXR observations from the gamma-ray spectrome- 
ter SPI onboard the Intern ational Gam ma-Ray Astrophysics 
Laboratory (INTEGRAL: iJean et al.llToOO,) were used. The 
energy bands of the RHESSI profiles were chosen in such a 
way that the emission was clearly nonthermal (steep rise - 
fast decay of the peaks). At the same time, the photon ener- 
gies had to be low enough to produce reasonable count statis- 
tics. Since the GOES classes of the analyzed events ranged 
between B and X, different energy bands of the RHESSI pro- 
files were used in each case to satisfy these demands. 

Data reduction: All images of a particular event were rotated 
to solar north, if necessary, and differentially rotated to the same 
reference time. Ha images were cross-correlated in time to ac- 
count for seeing effects. Coalignment of the different data sets 
was accomplished using MDI continuum. Ha blue wing, and 
TRACE WL images. Sunspots near the flare sites were taken 
as a reference for coalignment by cross-correlation techniques. 
The different pointings of the TRACE WL and 1600 A tele- 
scopes were taken into account. MDI images were converted 
from SOHO-view to Earth-view. 
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Fig. 1. 2005 January 15, X2.6 flare - a): Image frame corrupted by many particle hits, b): The same image after the particle hits 
have been removed by the running median method, c): Flare area at 22:46:43 UT detected from b). 



Table 1. Event Information. For each event we list the date, NOAA active region, position, and GOES class plus GOES 
maximum of the flare, as well as the first LASCO C2 appearance, the central position angle, and th e linear velocity of the 
associated CME. Details on the CMEs were obtained from the SOHO LASCO CME CATALOG dYashiro et all |2004 - 
http : //cdaw . gsf c . nasa . gov/CME JList/). 







Flare 






CME 




Date 


NOAA 


Position 




GOES 


first LASCO C2 


central PA 


linear velocity 




AR 




Class 


Max. [UT] 


appearance [UT] 


[deg] 


[kms-'] 


2003 Oct 28 


10486 


S16E08 


X17.2 


11:10 


10:54 


Halo 


2460 


2003 Nov 18 


10501 


S02 E37 


M3.9 


08:31 


08:50 


Halo 


1660 


2005 Jan 15 


10720 


N13 W04 


X2.6 


23:02 


23:06 


Halo 


2860 


2006 Jul 06 


10898 


S10W30 


M2.5 


08:35 


08:54 


Halo 


910 


2007 May 19 


10956 


NOl W05 


B9.5 


13:02 


13:24 


260 


960 



Table 2. Information on the used data sets. For each event we list the instrument, wavelength, cadence, and pixel scale of the image 
time series. A 'yes' in the Saturation-column indicates that the observations were saturated during the impulsive phase of the event. 
The last two columns give the used instrument and energy band of the observed nonthermal HXR emission. 



Date 




Imag 


; time series 






HXRs 


Instrument 


Wavelength 


Cadence [s] 


"/pixel 


Saturation 


Instrument 


Energy [keV] 


2003 Oct 28 


Meudon Obs. 


Ho- blue wing 


~ 60 


1.9 


no 


INTEGRAL/SPI 


> 150 


2003 Nov 18 


TRACE 


1600 A 


~ 23 


0.5 


no 


RHESSl 


20-60 


2005 Jan 15 


TRACE 


1600 A 


~ 10 


0.5 


no 


RHESSI 


30 -100 




BBSO 


Hq' line center 


~ 60 


1.05 


yes 




and 100 -300 


2006 Jul 06 


Hvar Obs. 


Ha line center 


~4 


0.3 


yes 


RHESSl 


20-50 


2007 May 19 


KSO 


Ha line center 


~ 60 


2.2 


yes 


RHESSI 


15-50 



When using intensity thresholds to differentiate between 
flare pixels and non-flare pixels, transient bright non-flare fea- 
tures, such as cosmic rays, will inevitably be included among 
ribbon pixels. This effect was distinct during the impulsive phase 
of one of the events studied (2005 Jan 15, X2.6). Since these 
features do not survive for more than 1 or 2 frames at any given 
pixel, they can easily be eliminated. One possibility is to cre- 
ate running mean images over a few frames. Thus, bright parti- 
cle hi ts are smoothed out (e.g.. iLongcope et ani2007l : lOiu et alj 
l2007h . although, the intensity of the flare pixels is also modified. 
Therefore, we created running median images in the image time 
series of this event. For each pixel x, we calculated the median 
of frames / to / + n (with n - 5) at the location of x. If x is a 
flare pixel, it will also be bright in the following frames and the 
median is high. If pixel jc is a particle hit, however, it will be dark 
again in the subsequent frames, and the median will be low. If 



the median is lower than a particular intensity threshold, pixel 
X in frame / will be replaced by the median, whereas all pix- 
els with medians exceeding the threshold remained unchanged. 
Thus, only particle hits were removed from the images, while 
the intensity of the flare pixels was not modified (cf. Fig. [Til. 
To check this method, we calculated the total flare area, derived 
from both the original time series of this event and the frames, 
where the particle hits had been removed, and found a total flare- 
area overestimation of ~100%, when using the original times se- 
ries. The overestimation of the reconnection flux, however, was 
only ~20%, since particle hits cover the entire FOV, i.e., also re- 
gions far away from the flare site, where the magnetic fields are 
weak. 
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3. Analysis 

To determine the magnetic-flux change rate ip - J B„da we 
measured the following quantities: 

1. The newly brightened area da in an image compared with 
the preceding images, separately for each magnetic polarity 
domain. This was done by: (1) creating base difference im- 
ages from the original time series; (2) determining the small- 
est intensity maximum Ism of the entire difference-image 
time-series; (3) multiplying Ism by various scaling factors 
to derive a set of potential intensity thresholds. Thresholds 
that enabled us to differentiate correctly between flare pix- 
els and non-flare pixels, were considered to be appropriate. 
The usability of the various threshold candidates was esti- 
mated by considering the total flare area detected in using 
them. Thresholds that were too low produced total flare areas 
that also contained non-flare pixels, whereas relatively high 
thresholds omitted fainter parts of the flare ribbons. Since 
several thresholds within a certain range seemed to be ap- 
propriate for identifying the total flare area, the analysis was 
carried out for a range of reasonable thresholds, and then the 
mean of the thus calculated magnetic-flux change-rate pro- 
files was taken. We found that the ^-peak-times were insen- 
sitive to the different thresholds, but the height of the peaks 
changed by 5 - 15%. 

To be counted as a member of da, a particular pixel had 
to exceed the given threshold, and be a non-flare pixel in 
the preceding images. Furthermore, it had to be located in- 
side the currently-analyzed magnetic -polarity domain, and 
exceed the MDI noise level of 20 G. 

2. The normal component of the photospheric magnetic-field 
strength B„ inside da. The line-of-sight magnetograms of 
MDI are known to be less sensitive to fields stronger than ap- 

Siroximat ely 1700 G, i.e., MDI underestimates strong fields 
B erger & Litesl2003h. Following the cross-calibration study 
of (Berger & LitesI (l2003h . we multiplied the reported MDI 
line-of-sight magnetic field values by 1.56. Furthermore, we 
assumed that the field is approximately radial at the photo- 
sphere, and therefore divided the photospheric magnetic field 
by the cosine of the central meridional distance of the flare 
to derive the radial magnetic-field strength at each pixel. 

At each time t, we measured both the newly brightened area 
da, consisting of all pixels meeting the aforementioned criteria, 
and B„ at the locations of these pixels, and then calculated the 
positive and negative magnetic reconnection flux at each time by 
multiplying the pixel area and B„ of each pixel and adding the 
products. Division of the reconnection flux by the time intervals 
between two consecutively taken images yielded the magnetic 
flux change-rate for the positive and negative polarity domain, 
(p+it) and (p-{t), respectively. The magnetic flux change-rate 
was calculated by taking the mean of <p+ and ip^. 

We also determined the cumulated magnetic reconnection 
flux at each time for both polarity domains, ip+{t) and ip-{t), re- 
spectively, by adding the newly reconnected flux at time f to the 
flux that had been reconnected up to time t. Since equal amounts 
of positive and negative magnetic flux are involved in the re- 
connection process at each time, the (p+{t) and (/3_(f)-profiles 
should be identical in the ideal case. We also calculated the to- 
tal reconnection-flux profile ^tot(0 to be the mean of (p+ and tp-, 
and we used the i/Jtot-profiles obtained from the lowest and high- 
est of the appropriate thresholds as an error estimate of the total 
reconnection flux. 



4. Results 

The results are presented separately for each event. The follow- 
ing is valid for all events. Panels (a) - (d) of Figs.|2]lll|6l[8l and 
[To] show the evolution of the flare ribbons in each event. Panels 
(a) were taken during the impulsive phase, panels (b) show the 
situation at the time of the maximum of the HXR profiles, pan- 
els (c) display the flare ribbons at the time of GOES maximum, 
and panels (d) are decay-phase images. Panels (e) show the cal- 
culated total flare area, i.e., the sum of the newly brightened ar- 
eas in all images, superposed on the decay-phase image, and in 
panels (f), the contours of this area are plotted on the MDI mag- 
netogram of the flaring region. The top panels of Figs. [3] |5] |7] 
E and [n] show the GOES 12 1 - 8 A soft X-ray (SXR) flux, as 
well as the cumulated reconnection-flux profiles, and in the bot- 
tom panels, the derived magnetic-flux change rate (p is presented 
along with the observed nonthermal HXR emission profile. 

4.1. 2003 October 28, X1 7.2 flare 

Panels (a) - (d) of Fig. |2]show the flare evolution in Meudon Ha, 
and panel (e) displays the total flare area. It consists of the two 
large regions, swept up by the separating flare ribbons, as well 
as some smaller regions in the vicinity that also brightened up in 
the course of the flare. In Fig.|2](panel f), the contours of the to- 
tal flare area are plotted on the MDI magnetogram. The northern 
ribbon passed through the negative magnetic polarity domain, 
whereas the southern one appeared in the positive polarity. Parts 
of both ribbons passed regions of magnetic field strength exceed- 
ing 1500 G. 



Meudon Observatory 28-Oct-2Cl03 11:03:14 Meudon Observatory 2B-Oct-2003 11:05:13 




-300 -200 -100 



100 -300 -200 -100 



Fig. 2. 2003 October 28, X17.2 flare - a) - d): Temporal evo- 
lution of the flare ribbons in Ha. a): impulsive phase image, 
b): time of INTEGRAL/SPI maximum, c): time of GOES maxi- 
mum, d): decay phase image, e): calculated total flare area super- 
posed on decay phase image (black with white contours: nega- 
tive polarity, gray with black contours: positive polarity), f): total 
flare area on MDI magnetogram. Contours are the same as in e). 
MDI data range scaled to ±1000 G out of [-1680, +1650] G. - 
FOV: 505" x 334". 
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11:00 11:05 11:10 11:15 11:20 11:25 

Start Time (28-Oct-03 10:57:00) 



Fig. 3. 2003 October 28, X17.2 flare - Top: GOES 12 1 - 8 A 
soft X-ray (SXR) profile, cumulated total magnetic reconnection 
flux (^tot) plus error estimate, and cumulated positive and nega- 
tive reconnection flux (^+, <^_). Bottom: INTEGRAL SPI HXR 
countrate and magnetic flux change rate {ip). 

The top panel of Fig.[3]shows the GOES flux, as well as the 
cumulated total reconnection flux along with its error estimate. 
For the sake of clarity, we plot 2 (/Jtot - ip+ + (f-.ln addition, the 
cumulated reconnection flux originating separately from the pos- 
itive and negative magnetic polarity domains, (/?+ and is pre- 
sented. The positive and negative flux profiles look similar over 
the entire time range. This indicates that almost equal amounts 
of positive and negative magnetic flux participate in the recon- 
nection process at a given time, as theoretically expected. During 
the impulsive phase, when more and more flux is reconnected, 
the cumulated flux profiles steeply rise. In the decay phase, i.e., 
after the GOES flux reached its maximum, the reconnection pro- 
cess slowly comes to an end, and the amount of newly recon- 
nected flux decreases. This results in nearly constant cumulated 
flux profiles during this phase of the event. At the end of the an- 
alyzed time interval, the ratio of cumulated positive versus neg- 
ative reconnection flux is 1.02, and the total flux (f^ot adds up to 
~ 15.5 X 10-' Mx (cf. Table©. 

The INTEGRAL SPI count-rate of > 150 keV, and the de- 
rived magnetic flux change-rate are presented in the bottom 
panel of Fig. [3] The INTEGRAL profile exhibits several spikes 
embedded in a broad peak, which lasts from ~ 11:02 UT to 
11:10 UT. Since the cadence of the Ha images was only about 
1 min, it was not possible to resolve every single INTEGRAL 
spike, although, the magnetic flux change-rate does reflect the 
overaU shape of the INTEGRAL flux. 

4.2. 2003 November 18, M3.9 flare 

This event was already published in iMiklenic et al. I (I2007h . 
We present it again here, because the cumulated magnetic 
reconnection-flux profiles had not then been calculated, but only 
the magnetic flux ratio of the overall positive to negative recon- 
nection flux, taken at the end of the analyzed time interval. 

Panels (a) - (d) of Fig. |4] show the flare evolution in 
TRACE 1600 A. In Fig. g] (panel e), the total flare area is 
presented, and in Fig. |4] (panel f) the contours of this area are 
superimposed on the MDI magnetogram. The northern ribbon, 
which sweeps across the negative polarity, stops when it reaches 
stronger fields. Thus, the flare area (white contours) is forced 
into a complex shape by the magnetic field topology. The south- 
ern ribbon is not affected by strong magnetic fields, so its shape 
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-100 
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-100 
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Fig. 4. 2003 November 18, M3.9 flare - Panels a) - d): Temporal 
evolution of the flare ribbons in TRACE 1600 A. a): impulsive 
phase image, b): time of RHESSI maximum, c): time of GOES 
maximum, d): decay phase image, e): calculated total flare area 
on decay phase image (black with white contours: negative po- 
larity, gray with black contours: positive polarity), f): total flare 
area on MDI magnetogram. Contours are the same as in e). MDI 
data range scaled to +500 G out of [-1700, +1220] G. - FOV: 
320" X 250". 

remains relatively simple, and the positive-polarity flare area 
(black contours) appears as an elongated structure. 

The cumulated positive and negative reconnection-flux pro- 
files look similar (cf. top panel of Fig.|5]l, although, slightly more 
negative than positive flux is detected near the end of the ana- 
lyzed time interval. The ratio of positive to negative flux is 0.91 
at about 08:30 UT (cf. Table |3]l. The total cumulated flux adds 
up to 2.3 X 10^' Mx. We note that the cumulated fluxes of the 
entire event may be higher than listed in Table [3] since it was 
not possible to analyze the event up to the end of the impulsive 
phase due to gaps in the TRACE data after ~ 08:24 UT. 

In the bottom panel of Fig.© the RHESSI HXR 20 - 60 keV 
time profile is presented along with the derived magnetic re- 
connection rate. During the event, RHESSI observed four major 
peaks. Three of them are also evident in the (^-profile, the time 
interval of a fourth peak being unable to be analyzed because of 
the aforementioned gaps in the TRACE data. The co-temporal 
occurrence of the first three RHESSI and ip-peaks is, however, 
obvious, the RHESSI peaks being slightly delayed, as indicated 
by the gray vertical bars, highlighting the time lag between the 
corresponding peak values. The numbers on the left of each bar 
indicate the time delay of the R HESSI peaks in seconds. We 
note that in ' Miklenic et al.l (12007 ). we did not divide the mag- 
netic field by the cosine of the central meridional distance of the 
flare, and therefore obtain in the present paper slightly higher 
values in the i/j-profile. 
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Fig. 5. 2003 November 18, M3.9 flare - Top: Same as in Fig.|3] 
The GOES-profile before 08:10 UT shows the occurrence of a 
previous flare (GOES class M3.2). Bottom: RHESSI HXR 20 - 
60 keV time profile and magnetic flux change rate (i^). Light- 
gray vertical bars mark time delays of RHESSI peaks compared 
to the associated magnetic flux change-rate peaks. Delay in sec- 
onds is given on the left of each bar. 



4.3. 2005 January 15, X2.6 flare 

Panels (a) - (d) of Fig. |6] show the flare evolution in 
TRACE 1600 A. Figure |6] (panel e) displays the total flare area, 
and, in Fig. |6] (panel f) the area contours are plotted superim- 
posed on the magnetogram of the flaring region. The negative- 
polarity area (white contours) is large compared to the positive 
one (black contours). This diff'erence in flare area is caused by 
the negative-polarity ribbon sweeping, at least partly, areas of 
weaker magnetic field, whereas the positive-polarity ribbon oc- 
cupies exclusively strong field regions, some of them exceeding 
1500 G. This produces comparable fluxes, although the size of 
the flare areas differ 

The horizontal bars in the top panel of Fig. |7] indicate 
the instrument observation times of BBSO and TRACE. Since 
TRACE observations did not start before ~ 22:42 UT, we used 
BBSO Ha images, which covered the entire impulsive phase, to 
estimate the cumulated reconnection fluxes for the time inter- 
val in which TRACE data was missing. We calculated the fluxes 
separately for Ha and TRACE 1600 A, and then we attached the 
flux profiles obtained from TRACE images to those derived from 
BBSO images up to the start of TRACE observations. Since the 
Ha images were saturated, the fluxes derived from BBSO may 
be overestimated, although, they still provide a reference level 
for the beginning of the TRACE profiles, and thus allow the es- 
timation of the overall cumulated flux in this event. The curve 
progression of the positive and negative fluxes is similar, yet, we 
detect more positive than negative flux during this event. At the 
end of the analyzed time interval, the ratio of cumulated posi- 
tive to negative flux is 1.35, and the total flux (^tot adds up to 
~ 4.5 X 10^' Mx (cf. TableO. 

The magnetic flux change-rate derived from TRACE images 
along with two RHESSI light curves are shown in the bottom 
panel of Fig. [T] The RHESSI peak between ~ 22:43 UT and 
22:53 UT is unusuafly broad, and in the 30 - 100 keV and 100 - 
300 keV profiles the highest RHESSI peak appears later than the 
maximum in (p. Until ~ 22:55 UT, the magnetic flux change rate 
exhibits four spikes embedded in a broad peak, while RHESSI 
exhibits the aforementioned broad main-peak. Thus, we did not 
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300 1 
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TRACE 1600 15^an-2005 22:43:50 4 TRACE 1600 15^an-2005 22:48:55 




TRACE 1600 15-Jan-2005 23:03:35 4 TRACE 1600 15-Jan-2005 23:09:34 




Fig. 6. 2005 January 15, X2.6 flare - Panels a) - d): Temporal 
evolution of the flare ribbons in TRACE 1600 A. a): impulsive 
phase image, b): time of RHESSI maximum, c): time of GOES 
maximum, d): decay phase image, e): calculated total flare area 
on decay phase image (black with white contours: negative po- 
larity, white with black contours: positive polarity), f): total flare 
area on MDl magnetogram. Contours are the same as in e). MDI 
data range scaled to +1000 G out of [-1470, -h1790] G. - FOV: 
160" X 134". 



assign the main RHESSI peak to any of the (p-peaks. The corre- 
spondence between the remaining RHESSI and (p-peaks is more 
obvious. After ~ 22:55 UT, the i^-profile and the progression 
of the nonthermal emission look similar, although the relative 
height of the RHESSI peaks is not always reproduced by the (p- 
profile, e.g., the last RHESSI peak at about 23:07 UT is high in 
comparison to its counterpart in the magnetic flux change rate. 
The co-temporal occurrence of observed and derived peaks is 
also obvious, however, in this event, the RHESSI peaks being 
slightly delayed. The vertical bars in the bottom panel of Fig. Q 
indicate the time lag between the peaks in the RHESSI 100 - 
300 keV and (^-profiles, and the numbers on the left of each bar 
indicate the time delay in seconds. We also checked the delay 
times of RHESSI peaks in the 30 - 100 keV profile and found 
delays of 53, 54, and 48 s. Both sets of delay times are compara- 
ble to those found in the 2003 November 18 event, presented in 
Sect.lO 



4.4. 2006 July 06, M2.5 flare 

Panels (a) - (d) of Fig. [8] show the flare evolution observed 
in Hvar Ha, which was associated with an erupting fila- 
ment. Figure [8] (panel e) displays the total flare area, and in 
Fig. [8] (panel f) the area contours are superimposed on the MDl 
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22:15 22:30 22:45 23:00 

Start Time (15-Jan-05 22:00:00) 

Fig. 7. 2005 January 15, X2.6 flare - Top: Gray and black hori- 
zontal bars mark the instrument observation times for BBSO Ha 
and TRACE 1600 A. Profiles derived from TRACE observations 
begin at ~ 22:42 UT and were attached to the light-gray profiles 
obtained from BBSO observations. Bottom: RHESSI HXR 30 - 
100 keV and 100 - 300 keV time profiles and magnetic flux 
change rate {(f). Light-gray vertical bars mark time delays of 
RHESSI 100 - 300 keV peaks compared to the associated mag- 
netic flux change-rate peaks. Delay in seconds is given on the 
left of each bar 



magnetogram. The negative -polarity area is small compared to 
the positive one, because the northern ribbon enters the penum- 
bra of the nearby sunspot, where the magnetic fields are strong. 
However, it stops at the border between penumbra and umbra. 

In this event, more positive than negative reconnection flux 
is detected over the entire analyzed time range (cf. top panel 
of Fig. |9]l. During the first part of the impulsive phase in par- 
ticular, the positive cumulated flux profile exhibits a steep rise, 
while the negative flux increases more gradually, its final flat- 
tening occurring later than that of the cumulated positive flux. 
Thus, cumulated positive and negative fluxes become compara- 
ble in strength at the end of the analyzed time interval, where 
the flux ratio is 1.06, and the total cumulated flux adds up to 
2.6 X 10^' Mx (cf. Table [3]l. We note that the Ha images were 
saturated during the impulsive phase of the event, so the calcu- 
lated flux values may be overestimated. 

In the bottom panel of Fig. |9l the RHESSI 20 - 50 keV 
time profile along with the derived reconnection rate are plot- 
ted. In this event, the nonthermal emission exhibits two dis- 
tinct peaks between 08:18 UT and 08:25 UT. Afterwards, sev- 
eral small spikes are visible between ~ 08:28 UT and 08:37 UT. 
The two main RHESSI peaks are not resolved in the magnetic- 
flux change rate. In the (/^-profile, only one broad peak is evi- 
dent, which spans the time range of the two main RHESSI peaks. 
The maximum of the broad (^-peak corresponds with the second 
RHESSI peak, and the RHESSI spikes occurring after 08:28 UT 
is evident in the magnetic-flux change rate as a broad bump. 
However, there is a second lower bump in the (^-profile after 
08:37 UT, where no corresponding increase in the nonthermal 
emission is visible. 

4.5. 2007 May 19, B9.5 flare 

Panels (a) - (d) of Fig. [TO] show the flare evolution observed 
in KSO Ha. There are two stable filaments east of the flare 
site. Two further filaments in this active region were erupting 
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Fig. 8. 2006 July 06, M2.5 flare - Panels a) - d): Temporal 
evolution of the flare ribbons in Ha. a): impulsive phase im- 
age, b): time of RHESSI maximum, c): time of GOES maxi- 
mum, d): decay phase image, e): calculated total flare area on 
decay phase image (black with white contours: negative polar- 
ity, white with black contours: positive polarity), f): total flare 
area on MDI magnetogram. Contours are the same as in e). MDI 
data range scaled to [-1000, -h400] G out of [-1390, -i-450] G. - 
FOV: 255" x 183". 




08:18 08:24 08:30 08:36 08:42 08:48 

Start Time (06-Jul-06 08:13:00) 



Fig. 9. 2006 July 06, M2.5 flare - Top: Same as in Fig.|3] Bottom: 
Magnetic flux change rate ((p) and RHESSI HXR 20 - 50 keV 
time profile. Counts at times of RHESSI shutter movements are 
set to 0. 

(cf. lVeronig et al.l2008l) . The first one, whose southern footpoint 
was anchored near the flare site, went off towards west and has 
akeady disappeared from the Ha filter at 12:50 UT. A part of 
the second erupting filament is still visible as a dark, elongated 
structure in panels (a) and (b) at the top of the images, although 
it is already in a state of eruption. In panel (c), it has almost left 
the Ha filter, and in Fig. [TO] (panel d) it has completely vanished. 
Figure[TO](panel e) shows the total flare area in this event, and in 
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Fig. [To] (panel f) the area contours are plotted on the MDI mag- 
netogram. The size of the positive and negative-polarity areas is 
comparable. Both areas are elongated structures. 
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Fig. 11. 2007 May 19, B9.5 flare - Top: Same as in Fig. [3] 
Bottom: RHESSI HXR 15 - 50 keV time profile and magnetic 
flux change rate {(p). 
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Fig. 10. 2007 May 19, B9.5 flare - Panels a) - d): Temporal 
evolution of the flare ribbons in Ha. a): impulsive phase image, 
b): time of RHESSI maximum, c): time of GOES maximum, 
d): decay phase image, e): calculated total flare area on decay 
phase image (black with white contours: negative polarity, gray 
with black contours: positive polarity), f): total flare area on MDI 
magnetogram. Contours are the same as in e). MDI data range 
scaled to +1000 G out of [-1900, -h1910] G. - FOV: 268" x 
306". 



In this comparatively weak event, more negative than posi- 
tive reconnection flux is detected (cf. top panel of Fig.fTTli. At the 
end of the analyzed time interval, the ratio of cumulated positive 
to negative flux is 0.64. The total flux adds up to 1.8 x 10^' Mx 
(cf. TableO. We note that the Ha images were saturated during 
the impulsive phase of the event, so the calculated flux values 
may be overestimated. 

The bottom panel of Fig. [TT] displays the derived magnetic- 
flux change rate, as well as the RHESSI 15 - 50 keV time pro- 
file. Although the lower end of this energy band could in prin- 



Table 3. Cumulated reconnection fluxes and ratio of positive to 
negative flux at the end of the analyzed time interval. Fluxes are 
given in units of 10^' Mx. Fluxes in the flare on 2003 November 
18, may be higher, since it was not possible to analyze the event 
up to the end of the impulsive phase due to gaps in the TRACE 
data after ~ 08:24 UT. Error estimates in i^tot are obtained from 
the lowest and highest of the appropriate intensity thresholds 
(see Sect.[3]for more detail). 



Date 


GOES 

class 




V- 




^+ 
If -I 


2003 Oct 28 


X17.2 


15.6 


-15.3 


15.5 


±0.8 


1.02 


2003 Nov 18 


M3.9 


> 2.2 


> -2.4 


> 2.3 


±0.2 


0.91 


2005 Jan 15 


X2.6 


5.2 


-3.9 


4.5 H 


=0.4 


1.35 


2006 Jul 06 


M2.5 


2.7 


-2.5 


2.6 H 


=0.4 


1.06 


2007 May 19 


B9.5 


1.4 


-2.2 


1.8 H 


=0.2 


0.64 



ciple be thermal emission, the shape of the RHESSI peaks and 
the RHESSI spectra indicate that the emission is predominantly 
nonthermal. The first and highest RHESSI peak is also evident in 
the (^-profile, whereas the second and third peak are unresolved. 

4.6. Relations between reconnection flux, GOES class, and 
CME velocity 

In Fig. [121 for the five flares studied we plot the measured GOES 
peak flux versus.: a) the total flare area A, i.e., the sum of the 
newly brightened areas in all images; b) the mean magnetic field 
strength (B) inside A; and c) the total magnetic reconnection 
flux i^tot- We calculated the three quantities for each of the ap- 
propriate intensity thresholds (see Sect. [3]l and plot the average. 
Error bars were obtained from A, (B ), and tptot for the lowest 
and highest thresholds. Error bars in Fig. [12] (panel b) also in- 
clude the effects of the MDI noise level of 20 G. After calcu- 
lating in log-log-space the correlation between the GOES peak 
of each event and A, {B ), and i^tot, respectively, the GOES peak 
was found to neither correlate with A nor {B ), although, it is cor- 
related significantly with i^tot (R = 0.89, confidence level 95%, 
cf. Fig.[T2](panel c)). Since our sample size is small, we note that 
we cannot exclude the possibility that A and (B ) are also related 
to the GOES peak flux (see Sect. l5.2l for a detailed discussion on 
this issue). 
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Fig. 12. GOES peak flux vs.: a) Total flare area A; b) Mean magnetic field strength (B ) inside A; c) Total reconnection flux. In the 
upper left corner of c) the Pearson correlation coefficient R, along with a, the probability of eiTor, is shown. EiTor bars are obtained 
from the lowest and highest of the appropriate intensity thresholds (see Sect. [3] for more detail). Error bars in (b) also account for 
the MDI noise level of 20 G. 



In Fig. [T3] (panel a), we combine the reconnection fluxes 
from the five events ana lyzed in this paper with r es ults from 
other events, derived by Qiu & Yur chvshvril (l2005h . lOiu et al] 
(|2007 |). and iLongcope etal. ( 2007) . We note that tiiese au- 
thors applied potential-field extrapolation from the photospheric 
MDI magnetograms to a height of 2 Mm, which we did not. 
Nevertheless, our results are consistent with the enlarged data 
set, and we derive comparable reconnection fluxes for the events 
on 2003 October 28 and 2003 November 18 (cf. encircled events 
in Fig.[T3Tl. The combined data set contains two events exceeding 
GOES class XIO (2003 October 28 and 29). For the entire data 
set, the correlation between GOES peak and i^tot is significant 
(R - 0.62, with a confidence level greater than 99%). However, 
if the two most energetic events are excluded, the correlation 
{R - 0.31) is significant only with a confidence level of 80%. 
This indicates that the rare > XlO-events strongly contribute to 
the correlation. 

We also investigated the relation between the total reconnec- 
tion flux in the flare and the kinematics of the associated CME. 
In Fig. [13] (panel b), we plot the lin ear CME velocity, ta ken from 
the SOHO LASCO CME catalog jYashiro et alj|2004 . i.e., the 
linear speed obtained by fitting a straight line to the height-time 
measurements, against i/Jtot for all events from Fig. [T3] (panel a). 
Using only the five events analyzed in this paper (marked as o in 
Fig. [T3] (panel b)), the correlation (R - 0.71) is significant only 
with an 80%-confidence level, although, taking the enlarged data 
set, we find a coiTelation that is significant with a confidence 
level greater than 99% (cf. Fig.[T3](panel b)), the coiTelation co- 
efficient (R - 0.76) being lower, however, than that reported by 
lOiu & Yurchvshviil (|2005|) for the subset of 13 events analyzed 
by these authors. Neverth eless, our result is in line with that of 
lOiu & Yurchvshvnl(l2005h . indicating that flares of higher recon- 
nection flux are associated with faster CMEs. 



03 V 
O ^ 

m 
m 
O 
C3 

M 



R = 0.62 (all events) 
a = 0.002 

R = 0.31 (allbu1X10.0,X17.2) 29oct-2003 i 
a = 0.20 



28-Oct-2003 



. o 

S 07-NOV-2Q04 



ie-Nov-2D03i;>a 

O 



a) 



R = 0.76 
a = 0.001 



28-Oct-2003 
lOAi 



ia-Nov-2003 i<>i' 



07-NOV-2004 



b) 



1 10 
Total reconnection flux [10^' Mx] 

Fig. 13. a) GOES peak flux vs. ^tot; b) Linear CME veloc- 
ity, taken from the SOHO LASCO CME catalog, vs. i^tot- <> 
mark event s from this study (cf. Fig. [12];): A indicate e vents 
taken from [Q iu & YurchYshYnl (l20()5l) r lOiu etaP i2QQf ). and 
[Longcope et a l. (2007). Same events analyzed in different stud- 
ies are encircled. Upper left corner of each panel shows Pearson 
correlation coefficient R along with a, the probability of error. 



5. Discussion 



5.1. Potential errors 



In the following, we list potential eiTors in the determination of 
the magnetic reconnection flux that may arise from the applied 
methods, as well as the technical equipment and data acquisition 
procedure. Afterwards, we discuss ffie physical implications of 
our results. 



Potentially significant measurement eiTors, which can result 
from transient bright non-ribbon features, such as cosmic rays, 
were accounted for by eliminating short-lived features from the 
images (see Sect. [2]for more detail). Disregarding transient fea- 
tures can result in an overestimation by 100% of the total flare 
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area and 20% of the total reconnection flux, if an image time se- 
ries is strongly corrupted by particle hits. Bright coronal loops 
may also overestimate the flare area, and thus, the reconnected 
flux, when counted among flare pixels. Therefore, we excluded 
these regions from the analysis by setting the corresponding pix- 
els during the time range, when these brightenings occurred, to 
zero. 

A crucial factor is the threshold value, which is u s ed to dis- 
cern flare pixels and non-flare pixels. lUongcope et al.l (l2007l) re- 
ported a change of 25% in the measured flux, when experiment- 
ing with cutoff values. We also tested several cutoff values and 
found a change of 5 - 15% in the measured reconnection flux. 

A further effect, which may cause an overestimation of the 
flare area, and thus, the reconnection flux, is caused by a limit to 
the amount of charge that each pixel of a CCD chip can store. 
If there is too much charge at a particular location on the CCD 
chip, i.e., if a pixel is saturated, it will overflow to its neighbor- 
ing pixels, preferentially the pixels above and below the satu- 
rated pixel. The signal in those pixels is questionable, since they 
will usually contain spilled charge. In three of the events (2005 
January 15, X2.6, 2006 July 06, M2.5, and 2007 May 19, B9.5), 
Ha images were saturated during the impulsive phase, and the 
area and reconnection flux may thus be overestimated. 

All sources of uncertainty mentioned so far are related to 
the flare area. The second quantity required to determine mag- 
netic reconnection rates and fluxes, namely, the magnetic field, is 
measured routinely only in the photosphere. Thus, photospheric 
line-of-sight magnetograms are commonly used for this purpose. 
However, since the magnetic field strength decreases from the 
photosphere to the chromosphere/transition region, where the 
flare area is measured, reconnection fluxes are overestimated, 
when the magnetic field is derived from photospher i c mag ne- 
tograms, as was done for the present paper lOiu et al.l (120071) re- 
ported an overestimation of about 20% in the reconnection flux 
derived from photospheric magnetic fields compared to the flux 
derived using potential field extrapolation to a height of 2 Mm. 

As for the magnetic fi eld, a further potential source of er- 
ror must be accounted for. iBerger & LitesI (l2003h reported that 
MDl measurements underestimate fields stronger than 1700 G, 
and Oiu et al. (2007) estimated an upper limit to possibly un- 
derestimated reconnection flux due to MDl saturation effects, 
assuming that saturation occurs at field strengths of 1700, 1600, 
and 1500 G. In two of the events analyzed in this paper (2003 
October 28, X17.2 and 2005 January 15, X2.6), small parts of 
the flare ribbons entered re gions with field strengths stronger 
than 1500 G. According to lOiu et all (l2007h . this may result in 
an underestimate of a few percent, at the very maximum, of the 
total reconnection flux. To ac count for this e f fect, w e followed 
the cross-calibration study of iBerger & LitesI ( l2003f) and multi- 
plied magnetic field values by a correction factor of 1.56. 

The misalignment between UV/Ho' images and MDl mag- 
netograms is another source of uncertainty, which can be es- 
timated by artificially offsetting the two sets of images. We 
did not experiment w ith i ntentionally misahgned i mages, but 
iLongcope etaTI (l2007h and lOiu & Yurchvshvnl (12005) reported 
that artificial misalignment by up to 4" can account for up to 10- 
20% of change in the measured flux. However, the co-alignment 
accuracy that we achieved was higher than 2". 

Bearing in mind all the sources of uncertainty discussed 
above, it is difficult to provide an overall error estimate for the 
total reconnection flux in each of the five events analyzed in 
this paper Effects of over- and underestimation may cancel each 
other to a certain degree. To be on the safe side, an overestimate 



of 30 - 40% should be taken into account, especially for events 
in which the fluxes were derived from saturated Ho' images. 

5.2. Consideration of tlie results in view of tlie standard 
flare/CME model and numerical MHD models 

We derived the magnetic-flux change rate from chromo- 
spheric/photospheric observations of five flare events (GOES 
classes B, M, and X) and compared it to the observed nonther- 
mal flare emission. In addition, we calculated the cumulated pos- 
itive and negative magnetic flux participating in the reconnection 
process, as well as the total reconnection flux iptot- We also deter- 
mined the correlation between the i^tot, GOES peak, and linear 
velocity of the flare-associated CMEs. 

We found good temporal correlations between the derived 
magnetic-flux change-rate and observed nonthermal emission in 
all events, i.e., hard X-ray peaks were clearly reflected in the 
magnetic flux change-rate profiles, as expected from the stan- 
dard model, although the relative height of the peaks was not 
always reproduced in the i^-profiles. In two events, for which 
non-saturated TRACE image time series were available with par- 
ticularly high cadence (< 23 s), the HXR peaks appeared de- 
la yed by roughly 1 min . Considering the amount of this delay, 
in iMiklenic et al.l (l2007h . we speculated that it might be related 
to the travel time of a reconnected field line from the diffiision 
region to the lower edge of the current sheet. However, since we 
found a delay of HXR peaks in only two events so far, no con- 
clusions can be d rawn a t prese nt. It is, however, interesting to 
note that Warmuth et al.l ( l2009l) calculated a travel time of 92 s 
in the 2003 November 18 flare, when comparing their model of 
shock drift acceleration at the reconnection outflow termination 
shock with observations. 

The derived total reconnection fluxes ranged between 1.8 x 
10^' Mx for the weakest event and 15.5 x 10"' Mx for the most 
energetic one (cf. Table [3]l. According to the standard model, 
equal amounts of positive and negative magnetic flux should 
participate in the reconnection process. The cumulated posi- 
tive and negative fluxes that we found were roughly balanced, 
with flux ratios ranging b etween 0.64 and 1.35. According to 
iOiu & Yurchvshvnl (120051) . flux ratios of between 0.5 and 2 can 
be regarded as a good flux balance, given the numerous uncer- 
tainties involved in the measurements. 

Beside determining reconnection rates and fluxes, we also 
calculated the correlation between GOES peak flux and: (1) to- 
tal flare area; (2) mean magnetic field strength inside this area; 
and (3) total reconnection flux. Although neither the flare area 
nor the magnetic field were correlated with the GOES class, the 
total reconnection flux was. However, we note that our sam- 
ple size is small. Therefore, we cannot exclude the possibil- 
ity that the other two parameters are also related to the GOES 
peak flux, and statistical studies have indeed found such corre- 
lations. Nagashima & Yokoyama (2006) analyzed 77 C, M, and 
X flares and detected a threshold value in the magnetic field that 
increases with the GOES peak flux. For example, X-class flares 
occurred only when the mean ma gnetic fie l d in th e active re- 
gion was stronger than ~ 100 G. ISu et aP (l2007l) analyzed a 
sample of 3 1 flares and found that for events of stronger aver- 
age magnetic field strength, the GOES peak flux tended to be 
higher. In a subsample of 18 out of the 31 events, they also de- 
tected a positive correlation between the GOES class and both 
the flare area and the magnetic flux, the magnetic flux show- 
ing a much st r onger correlation than the area or the magnetic 
field. ISu et al.l (l2007l) attributed this to the magnetic flux being 
the product of the other two parameters. Although the applied 
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methods for determining the flar e area differ from our approach - 
iNag ashima & Yokovama' ( 2006) used intensity thresholds in soft 
X-ray images, and Su et al. ( 2007) took a contour level in photo- 
spheric magnetograms of the flaring region as a basis for the cal- 
cu lation of A - our findings are in line with the explanation given 
bv lSu et al] (|2007|) . because they also indicate that in this context 
the combination of area and magnetic field inside this area is im- 
portant. It is conceivable that we did not find the weaker corre- 
lations between GOES class and A and {B ) due to the small size 
of our sample, but we did find the stronger correlation between 
GOES peak flux and (^tot- 

In addition, we incorporated our flux results in t o a larger 
set of result s, taken from Oiu & Yurchvshvnl (|2005|) . lOiu et al] 
(|2007|) . and iLongcope et alT ( 20071) . and again found a signifi- 
cant correlation between GOES peak flux and total reconnection 
flux (R - 0.6, confidence level greater than 99%). However, the 
correlation became smaller (R = 0.3), and was significant only 
with an 80%-confidence level, when the most energetic events 
(GOES classes > XIO) were excluded. This suggests that the 
rare, most energetic flares strongly contribute to the correlation. 
Still, the results indicate that the amount of magnetic flux par- 
ticipating in the reconnection process is larger in more energetic 
events than in weaker ones, as theoretically expected. The more 
magnetic flux is reconnected in a flare, the more energy is re- 
leased into fast particles and can subsequently be deposited into 
the chromosphere. This energy heats the chromospheric plasma, 
which then evaporates, fills the flare loops, and causes them to 
emit soft X-ray radiation, measured by the GOES satellites. The 
event is assigned to a particular class in the GOES classification 
scheme according to its maximum soft X-ray emission, which 
can be understood as a measure of the cumulated e nergy in the 
hot thermal plasma (e.g.. lVeronig et al.ll2002ll2005l) . 

We also found that flares with more reconnection flux are 
associated wit h faster CMEs. This co rrelation, which was also 
reported by Oiu & Yu rchvshvril ( 12005 ). can be explained by a 
feedback relationship between the CME kinematics and the re- 
connection process in the fl are (IVrsnak et al.l2005l : lMaricic et al.l 
120071; iTemmer et al.l200"8l) . According to the standard flare/CME 
model, 50% of the magnetic flux reconnected during a dynam- 
ical (two-ribbon) flare is transported upward to higher coronal 
heights and even interplanetary space. The reconnection process 
in the flare, which is initiated in the wake of the CME, affects the 
CME by: (1) reducing the net tension of the overlying magnetic 
field; (2) increasing the magnetic pressure below the erupting 
flux- rope; and (3) supp lying the flux-rope with extra poloidal 
flux jMaricic et al J 12007) . This enhances and prolongs the ac- 
celeration of the flux-rope and therewith additionally drives 
the CME. Therefore, CMEs that are associated with dynamical 
flares are generally faster than those linked with eruptive fila- 
ments, and CMEs that are associated with more powerful flares, 
i.e., flares of higher GOES class orfluence (time-integrated SXR 
flux), are of higher medi an speed and kine t ic energy than weak - 
flare-associated CMEs dMoon et all 120021; IVrsnak et all llOOl) . 
In other words, the higher the amount of flux reconnected in a 
flare, the higher the amount of flux going upward, resulting ulti- 
mately in faster CMEs. 

Numerical MHD models deaUng with the formation of a flux 
rope, its sudden eruption and escape from the Sun provide ad- 
ditional evidence of a feedback relation ship between an er upt- 
ing structure and magnetic reconnection. lCheng et alj (l2003h ac- 
complished 2.5-dimensional resistive MHD simulations of the 
evolution of an initially closed arcade field configuration. By 
imposing a shear-increasing footpoint motion, magnetic recon- 
nection occurs, creating a new flux rope. The toroidal flux origi- 



nally contained in the arcade is redistributed into both the slowly 
rising flux rope and the underlying arcade field. Thus, the mag- 
netic shear in the arcade is reduced after the flux rope formation, 
but increases again, if additional footpoint motion occurs. As 
soon as the magnetic shear exceeds a critical value, another re- 
connection event occurs and a new flux rope is formed. When 
magnetic reconnection continues, this newborn flux rope rises 
with increasing velocity, and finally merges with the previously 
created flux rope. This process of flux rope formation and merg- 
ing can be repeated as long as magnetic shear is continuously 
supplied. ICheng et alj (|2003[) quantitatively compared the model 
results with low-corona observations of a filament eruption dur- 
ing a flare/CME event and found good agreement in the evolu- 
tion of the flux rope height, velocity, and acceleration during the 
flux rope acceleration phase. Although the model does not ad- 
dress the escape of a flux rope from the Sun into interplanetary 
space, it shows that reconnection enhances the flux rope acceler- 
ation, and thus supports the scenario of a feedback relationship 
bet ween the erupting str ucture and the reconnection in the flare. 

iLin&ForbesI (120001) investigated the interaction between an 
already existing flux rope and magnetic reconnection. In their 
two-dimensional model, the flux rope experiences a catastrophic 
loss of equilibrium when the photospheric sources of the coro- 
nal magnetic field are brought together quasi-statically and reach 
a critical distance. Then, the flux rope jumps to a new equilib- 
rium position at a higher altitude, and a vertical current sheet 
is created below the flux rope. The model shows that without 
reconnection in the current sheet the magnetic tension force is 
always strong enough to prevent the flux rope from reaching in- 
terplanetary space. However, even a modest reconnection rate 
(prescribed by the inflow Alfven Mach number Ma) is sufficient 
to allow its escape. Based on this loss-of-equilibrium model and 
its exten sion, whic h includes g ravity i n the calculations dLinl 
2004; Re eves & Fo rbes 2005b), Reevesl ( l2006l) examined the rel- 
evance of reconnection for the flux-rope acceleration by analyz- 
ing its impact on the individual forces acting on the flux rope. 
They found that only the force due to the current sheet was 
affected considerably by changes in the reconnection rate Ma, 
while other forces, e.g., gravity, exhibited only minor changes 
with varying Ma. Slow reconnection rates resulted in longer cur- 
rent sheets, and thus, an increase in the downward force exerted 
by the current sheet on the flux rope. As a consequence, the total 
acceleration of the flux rope was decreased in cases of slow re- 
connection rates. When the reconnection rate was fast, the force 
due to the current sheet was diminished, because fast reconnec- 
tion rates dissipate the current sheet. Therefore, the total flux 
rope acceleration was higher in cases of fast reconnection. 

All of the models aforementioned are 2D or 2.5D models, 
which might be appropriate only in describing the eruption of a 
very long flux rope, especially when its height and width are 
much smaller than its length. In 3D models, the flux rope is 
anchored in the photosphere, and because of its final size, the 
overlying field can be pushed aside by the erupting structure. 
Thus, unlike in the 2D-approach, in 3D-model reconnection is 
not necessarily the main mechani s m red ucing the tension of the 
overlying field. However. I VrsnakI (l2008h demonstrated that in a 
3D-flux-rope model without reconnection high accelerations are 
prevented because of the inductive decay of the current in the 
flux rope as the eruption progresses. On the other hand, includ- 
ing reconnection underneath the rising flux rope provides 'fresh' 
poloidal flux to the rope. This preserves the current in the flux 
rope, and thus prolongs and enhances its acceleration. 

The theoretical work discussed above reveals the importance 
of magnetic reconnection for the kinematics of the erupting 
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Structure. Reconnection in the current sheet occurs during the 
eruption of the flux-rope/CME, additionally drives the eruption, 
and may even be essential for a flux-rope to escape from the Sun, 
i.e., there may be no fast CMEs without magnetic reconnection, 
which is generally believed to play a major role in solar flares. 
However, it is known from observations that CMEs are not nec- 
essarily associated with flares. Instead, they may occur together 
with erupting filaments. After the eruption of a quiescent fila- 
ment a growing system o f post-eruption loop s is often observed 
in the EUV range (e.g.. IVrsnak et al]|2005h . Morphologically, 
these loop systems are similar to postflare loops, except for 
the fact that they are not hot enough to be seen in soft X- 
rays. Quiescent filament eruptions occur in quiet regions, where 
the plasma-to-mag netic-pressure ra ti o (i is generally larger than 
in active regions dWu et all llOOlh . IVrsnak & Skendej (l2005h 
demonstrated that in a yS > 0.1 environment, the reconnection 
outflow region is not heated much. Therefore, an eruption taking 
place in such an environment should not be accompanied by a 
flare. 

The loss-of-equilibrium model also indicates that insuf- 
ficient plasma heating r nay be responsible for non-flare 
CMEs. iReeves & ForbesI (l2005ah simulated several erupting 
flux-rope/CME events, using background magnetic fields that 
ranged from weak to higher fields strengths. They found 
that weaker/higher background fields resulted in slower/faster 
CMEs, if all other parameters in the calculations remained the 
same. They also calculated GOES X-ray light curves for the 
1 - 8 A channel adding two diff'erent levels of soft X-ray back- 
ground flux, namely, the average of the maximum/minimum 1 - 
8 A GOES flux during solar cycle 22. It turned out that with 
the weakest of the chosen background magnetic fields, adding 
the low X-ray background level produced an A5 flare, while 
adding the high X-ray background-flux, the GOES light curve 
hardly emerged from the background, i.e., in this case, the ther- 
mal energy release was not large enough to heat and evaporate 
sufficient plasma into the loops. Thus, the corresponding GOES 
light curve progression should probably not be considered to be a 
flare, although, magnetic reconnection was involved in the erup- 
tion of the flux rope/CME. 

It is also known from observations that fa st CMEs are 
not a lways associated with powerful flares (e.g, IVrsnak et alj 
12005*). Although there are correlations between the kinetic en- 
ergy/velocity of CM Es and the SXR peak flux in the corre - 
sponding flares (e.g.. iMoon et alj|2003l: [Surkepile et al.ll2004l) . 
CMEs of comparable speed can by all means be associated with 
flares of different GOES classes. Two of the flare/CME events 
analyzed in this paper provide an observational example. The 
CME-speeds on 2006 July 06 and 2007 May 19 were of the or- 
der of 900km s ', but the GOES importance of the associated 
flares was M2.5 and B9.5, respectivel y (cf. Table [T]i. The sim- 
ulations carried out bv ' Reeves & Forbe s (2005a) also addressed 
this issue, and the results indicated that the CME speed is not 
indicative of the amount of associated X-ray flux, i.e., CMEs 
with very similar trajectories can have quite diff'erent flare re- 
sponses. The authors simulated a high-mass flux rope erupt- 
ing from a strong background field region as well as a low- 
mass flux-rope being ejected out of a weaker background field. 
The resulting velocity profiles for both events were similar The 
flux-rope/CME speeds reached values of around 1100 km s"', al- 
though the corresponding GOES peak values, and thus, the im- 
portance of the associated simulated flares were diff'erent (C5 
vs. M2).|R eeves & ForbesI (l2005al) attributed the difference in 
GOES class to the amount of thermal energy released in the cur- 



rent sheet in the model being directly related to the strength of 
the background magnetic field, while the kinetic energy of the 
flux rope depends on both the mass of the flux rope and the 
strength of the background field. Therefore, low mass CMEs 
from weak-field regions can have a similar velocity profile as 
massive CMEs from stronger field regions. The difference in the 
physical properties of both situations is evident only in the sim- 
ulated X-ray emission of the associated flare. 

Taking into account the results provided by numerical MHD 
simulations of solar eruptions, it becomes apparent that the com- 
bination of theoretical work and flare/CME observations is es- 
sential for improving our understanding of the underlying phys- 
ical mechanisms and the relation between the two phenomena. 
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